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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-1666 

A GENERAL DIGITAL COMPUTER ANALYSIS OF 
STATICALLY INDETERMINATE STRUCTURES 1 
By Paul H. Denke 

SUMMARY 


The application of high speed digital computers in the rational analysis of 
statically indeterminate structures, and the significance of this application in 
airframe design, are discussed. 

The matrix • formulation of the force method of analysis is reviewed, and the 
programs which have been produced to generate the matrices and solve the equilib- 
rium and continuity equations are described. These programs are general enough 
to apply to any linear discrete structure. 

Numerous comparisons between analysis and experimental results are presented. 
In addition, applications of the programs in the production stress analysis of a 
large commercial jet transport are described. Applications to thermal stress 
problems and low aspect ratio wings are also included. 


llhis paper, which carried a Douglas Aircraft Company designation of 
"Engineering Paper No . 83 L, " was presented before a meeting of the Structures and 
Materials Panel of the Advisory Group for Aeronautical Research and Development, 
North Atlantic Treaty Organization, in Aachen, Germany, September 17, 1959- Since 
the proceedings of the above Panel meeting are not being published, arrangements 
have been made with AGARD and the Douglas Aircraft Company for the release of this 
paper in its original form by NASA to increase its availability. 



NOTATION 


In the following definitions, the terra "analysis condition" means any com- 
bination of external load, thermal deformation, support displacement, etc., tend- 
ing to produce stress and deflection in the structure. The matrices are defined 
in the order of their appearance in the analysis. Matrices which are not in the 
liBt are defined in the teKt. 


Matrix 

• -M 


1 -M 
■ M 

' -M 

% ■ [”pp 13 ] 


Definition of the Matrix Element 
* the ith principal statically determinate force result- 
ing from the redundants and the external loads in the Jth 
analysis condition. 

Q - the ith subordinate statically determinate force re- 

8 iJ 

suiting from the redundants and the external loads in the Jth 
analysis condition. 

Xij » the ith principal redundant in the Jth analysis condi- 
tion. 

X ■ the ith subordinate redundant in the Jth analysis con- 

8 iJ 

dition. 

m the ith principal external load in the Jth analysis 
condition. 

0 ■ the ith subordinate external load in the Jth analysis 

8 iJ 

condition. 

m • the component in the ith principal degree of freedom 

of a unit value of the Jth principal statically determinate 


force . 
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Definition of the Matrix Element 



m = the component in the ith principal degree of freedom 
pS iJ 

of a unit value of the Jth subordinate statically determin- 
ate force. 

p ss -che component in the ith principal degree of freedom 

* pp iJ 

of a unit value of the Jth principal redundant. 


the component in the ith principal degree of freedom 


PS 


U 


of a unit value of the Jth subordinate redundant. 


the component in the ith principal degree of freedom 


PP 


ij 


of a unit value of the Jth principal external load. 


the component in the ith principal degree of freedom 


ps 


U 


of a unit value of the Jth subordinate external load. 


m 


sp 


» the component in the ith subordinate statically de~ 


ij 

terminate degree of freedom of a unit value of the Jth prin- 
cipal statically determinate force. 

m = the component in the ith subordinate statically de- 
BS iJ 

terminate degree of freedom of a unit value of the Jth sub- 
ordinate statically determinate force. 


the component in the ith subordinate redundant de- 


ep. 


U 


gree of freedom of a unit value of the Jth principal redun- 
dant. 


the component in the ith subordinate redundant de- 


ss 


ij 


gree of freedom of a unit value of the Jth subordinate redun- 
dant . 


the component in the ith subordinate external load 


ep 


U 


degree of freedom of a unit value of the jth principal exter- 
nal load. 
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Matrix 


Definition of the Matrix ELement 



p - the component in the ith subordinate external load 
8B iJ 

degree of freedom of a unit value of the Jth subordinate 
external load. 


k 



INTRODUCTION 


For many years, elementary methods of Btress analysis were used almost ex- 
clusively in the design of aircraft structures. These methods involved a number 
of assumptions , including especially the assumptions that plane sections of e* 
longated members remained plane under the action of bending loads, and that, in 
torque, sections were free to warp. In many parts of the airframe these assump- 
tions were, and are, completely Justified by the nature of the structure and the 
loading. In other places, the assumptions did not apply, as at the roots of 
wings, or in the regions of fuselage cutouts. In such areas, other assumptions, 
conservative and often overlapping to ensure safety, were made. Occasionally a 
more precise analysis was performed, but such occasions were rare. 

Actually no other recourse was possible, because the extensive use of pre- 
cise methods required computing facilities which did not exist. Such facilities, 
however, are now available. To appreciate the advance which has been made in 
the art of computation, consider the fact that about twenty seconds are required 
to multiply two seven digit numbers on a desk calculator, whereas a large auto- 
matic computer can multiply 10,000 pairs of such numbers per second. These fi- 
gures represent an increase in computing power on the order of 200,000 : 1. On 
a cost basis, the expense of computing has decreased on the order of 5,000 : 1. 

The introduction of matrix algebra into structural analysis has facilitated 
calculations also, by converting what was formerly a complicated mathmatical pro- 
blem into a systematic procedure. 

The result of these improvements 1 b that the use of advanced methods in 
stress analysis is now a practical undertaking. The question is, to what extent 
should these methods be applied. 

Figure 1 shows the results of a test run at NASA on a cylindrical shell sup- 
ported at , one end on a rigid foundation, reinforced by circular rings, and carry- 
ing a radial load at the free end. The figure shows the longitudinal tensile and 
compressive stresses in the shell, as determined from test, as computed by ele- 
mentary theory (My/l), and as computed by rigorous methods. The figure shows that 
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the maximum bending stress at station 45 frame as computed by elementary theory 
is in error by a ratio of almost 3*6 to 1, whereas the error resulting from the 
rigorous computation is only 10/., Notice also that a secondary maximum occurs 
at the so called "neutral axis" where the stress is supposed to be zero. Even 
at the rigid support, where the section is forced to remain plane, the error in 
My/l is still 2.2 to 1. This structure is not an Isolated case; it is typical 
of many parts of the airframe, and there are places in actual structure where 
errors resulting from elementary analysis may be larger, because of the exis- 
tence of cutouts or other conditions. 

The results of Figure 1 are well confirmed, inasmuch as they were obtained 
independently by Jensen of the Gruman Aircraft Company and published by him in 
reference 5 * These resultB cannot be Ignored or dismissed} they are facts, and 
must be considered in any assessment of structural analysis methods. 

What is the significance of the errors Involved in the use of elementary 
methods? 

Structure analyzed by rough methods and not thoroughly checked by a care- 
ful testing program can contain large stress concentrations. These concentra- 
tions can produce metal fatigue and cause the structure to have a short life. 

Much importance has been attached, justifiably, to the effects of small scale 
stress concentrations around bolt holes, tool marks, small radius fillets), etc., 
in reducing fatigue life. Perhaps not enough emphasis has been given to the im- 
portance of large scale stress concentrations that are not revealed by rough 
analysis methods. Obviously, an unconservative error of 3 t 1 or more in the 
computed stress, if undetected, must lead to a short lived structure. In such 
a case no amount of attention to design details, important as they are, can pro- 
duce a fatigue resistant component. The possibility exists that many of the fati- 
gue troubles experienced in the operation of present day aircraft have resulted 
from the use of elementary stress analysis methods where they did not apply. 

These large scale stress concentrations can also cause failure under the 
action of a single load, even though yielding tends to alleviate the condition. 

The consequences of such a failure need not be emphasized. 

If, as is normally the case, a thorough testing program is undertaken, then 
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all stress concentrations of importance can be discovered and eliminated. How- 
ever the cost of building, instrumenting, and testing full scale components is 
very high, even compared to the rental of a large computer. This testing ex- 
pense continually increases as the demand for higher performance vehicles re- 
quires the working of metals to higher operating stresses, the use of unusual 
configurations, and the ability to withstand severe environmental conditions. 

The testing of large components and entire airframes at high temperature will be 
an especially expensive procedure, because of the large power requirements to 
heat, as well as to cool, the specimen; the complicated apparatus needed for tem- 
perature control; the specialized instrumentation, such as high temperature strain 
gauges required for measurements; and the additional engineering required to plan 
the test. The new methods of stress analysis can play a very important part in 
helping to keep these testing expenditures within reasonable limits. 

Finally, the financial risk involved in a large aircraft project is suffi- 
cient to warrant a double check through both test and accurate analysis to make 
sure that no defective conditions exist. 

The conclusion is drawn, therefore, that the extensive use of advanced digi- 
tal methods of stress analysis is Justified at the present time, and that these 
methods will become even more important in the future. 

SCOPE OF THE PAPER 

The paper contains a general description of the method and sections on the 
matrix formulation, computer programs, analysis procedures, comparisons with test 
results, and applications. For a non-technical description of the work, the sec- 
tions on the method, test results, and applications are recommended. 
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TIIE METHOD 

In the following discussion, the terra "discrete structure" denotes a struc? 
ture composed of a finite number of members connected at a finite number of 
joints. The term "linear structure" denotes a structure for which the relation 
ships between external load, support displacement, internal force, and deflec- 
tion are linear. 

Almost every procedure for the analysis of statically indeterminate struc- 
tures can be classified as either a "force" or a "displacement" method. In the 
force method, the unknown internal forces are calculated first; the displace- 
ments second. In the displacement method, the displacements are calculated be- 
fore the forces. Argyris [l] * has discussed the two methods and shown the exist- 
ence of an analogy between them. 

The capabilities of the digital computer allow either of the basic method* 
to be programmed in ita simplest and most general form. In the past, a great 
many variations of the basic methods have been employed. One reason for such 
diversity has been the need to avoid extensive calculation by tailoring the 
method to fit the structure. However, the development of the digital computer 
has altered the situation. Extensive calculations now can be performed rapidly 
and economically. Therefore, a return to basic principles is feasible and, 
furthermore, the computer program designed to utilize these principles can be 
general in its applications . 

Some of the advantages to be gained from a basic, general approach are re- 
duced programming time, reduced training of personnel, the added insight that 
results from the application of basic principles, and the reduction of errors 
that results from familiarization in the use of a single method. 

The method of analysis described in this paper is a matrix formulation of 
the equilibrium equations and the Maxwell -Mohr equations for statically indeter- 
minate structures. This formulation was presented at a meeting of the Second 
U.S. Congress of Applied Mechanics in June, 195^ [XI • The use of matrix algebra 
is now recognized as essential in preparing the structural analysis problem for 
the computer. Langefors £ 3 ] and Wehle and Lansing [VJ had previously published 

* Numerals in brackets indicate references. 
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matrix formulations of Castigliano ' s Theorem. However, the Maxwell -Mohr equa- 
tions are a little simpler in form because they do not involve partial deriva- 
tives. Also, the applications to thermal stress and nonlinear problems are 
more straight-forward. 

In the Maxwell -Mohr method, which is a force method, the structure is cut 
to create a statically determinate structure or basic system. The members of 
the statically determinate structure may be simple elements, or they may them- 
selves be complicated statically indeterminate structures, (in fact, even so 
called simple elements are actually Infinitely redundant). After cutting, values 
of the redundants are chosen such that the deflections at the cuts resulting 
from external loads, support displacements, element thermal and other deforma- 
tions, and from the redundants, are zero. The redundants can be either forces 
existing at the cuts, or linearly independent combinations of these forces, as 
Argyris has pointed out [Y]. The conditioning of the simultaneous equations in- 
volved in solving for the redundants can be improved either by cutting on the 
basis of physical reasoning so that the forces at the cuts are small compared 
to other forces in the structure, or by linearly transforming the redundants 
on the basis of the known orthogonal solution of a geometrically regular struc- 
ture which bears a resemblance to the structure under consideration. The use of 
statically indeterminate substructures as elements, which have been previously 
analyzed, also improves the conditioning. 

The present method comprising the equilibrium and Maxwell -Mohr equations 
and the associated digital computer program is applicable to any linear discrete 
structure, and through iterative techniques to certain nonlinear structures as 
well. The method applies not only to various parts of the airframe structure 
such as the wing -fuselage intersection, the tail-fuselage intersection, the cock- 
pit enclosure, the area surrounding a fuselage cut-out, a low aspect ratio wing, 
and so on, but also to many types of structures encountered in civil engineering 
practice . 

This generality was not designed into the method to show the versatility 
of the computer, but because generality is necessary if the analyst is to have 
the tools that he needs to deal with the problems arising in airframe and missile 
design. Thus, many important airframe components have no recognizable geometric 
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regularity such as would permit the use of simplifying but restrictive assump- 
tions, or the application of results from elasticity theory. Figures 2a and 2b, 
which show a pylon -wing intersection, illustrate a structure of this kind. 

MATRIX FORMULATION 

The matrix formulation is preceded by a set of equations in vector notation 
which permit the calculation of the elements of the equilibrium matrices. 

Eq uil ibrium equations for a statically determinate structure are written by 
setting the sum of components of forces in a given direction and the sun of mo- 
ments about a given axis equal to zero. In general, such a set of equations can 
be expressed in matrix notation in the form MQ + P 0 <= 0. In this equation, Q 
is a matrix of unknown generalized forces where the term "generalized force" iB 
understood to mean either a force or a moment. The coefficients of the unknown 
forces Q are contained in M. These coefficients, called generalized components, 
are force or moment components in certain directions or about certain axes of 
unit values of the generalized forces. 

The matrix is a matrix of external loads acting on the structure, while P 
contains generalized components of unit values of these external loads. 

The structure to be analyzed is broken into free bodies, and equilibrium 
equations are written for each body. The equations are numbered consecutively 
beginning with one, and to each equilibrium equation there is assigned a corres- 
pondingly numbered unit vector coinciding with the direction in which forces are 
summed or about which moments are taken. These vectors are called degree of free- 
dom vectors , because only as many of them may be assigned to a free body as the 
body has degrees of freedom if the corresponding equations are to be independent. 
Figure 3 shows a free body diagram with forces and degree of freedom vectors re- 
presenting equations of equilibrium. Degree of freedom vectors are shown dotted. 

The existence of two kinds of equilibrim equations and two kinds of general- 
ized forces means that there can be four kinds of generalized components. Equa- 
tions 1, 2, 3, and 4 of Table 1 provide the method for calculating these quan- 
tities. In these equations, T^ is a unit degree of freedom vector (either trans- 
lational or rotational), and Fj Is a unit generalized force (either a force or a 
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moment). The symbol ra^j denotes the corresponding generalised component . In the 
rotation-force equation, r^ is a vector Joining the origin to any point on the 

lino of action of Tj, and fj is a similar vector Joining the origin to any point 
on the line of action of Fj. In equations (l) to (4), the frame of reference le 
assumed to be a right-handed rectangular Cartesian coordinate system, and rota- 
tions and moments are represented by vectors according to the right-hand rule. 

After the statically indeterminate structure is cut, three kinds of forces 
are seen to be acting upon, or in, the determinate structure. These forces are 
the external loads, the redundants, and the unknown internal forces, referred to 


TABLE 1 

SUMMARY OF EQUATIONS 


GENERALIZED FORCE COMPONENTS 


Translation-force 

m iJ 

" T i 


(1) 

Rotation-force 

m iJ 

- T i 

• [( ; J- r i ) * F J 

(2) 

Translation-moment 

m iJ 

ts 0 


(3) 

Rotation-moment 

“ij 

- T i 


(4) 


THE K TRANSFORMATION MATRICES 




-l 

■ -in m 

(5) 

ss sp 

-1 

m — p p 

X X„ 

(6) 

ss sp 


-1 

*» -p p 

*o 

(7) 


K - 


ss sp 

COEFFICIENT MATRICES IN THE PRINCIPAL EQUILIBRIUM EQUATION 


M 


m + m K 
pp ps m 


(8) 
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(9) 


P ■ p + p K 
x *x *x a 
PP PS 


P. * P rt ♦ P. K 

o o o o 

PP PB 


( 10 ) 


STATICALLY DETERMINATE FORCES RESULTING FROM UNIT REDUIDiWm AND UNIT EXTERNAL 
LOADS 


- -M" 1 P x ( 11 ) 

% “ -M' 1 P q (12) 

ELEMENT FORCE AND STATICALLY DETERMINATE REACTION MATRICES 


f x - 

N «x + H X 

(13) 

r Dx 

" N r«x 

(16) 

f - 

o 

N «o + H o 

(14) 

r Do 

" N r<o 

(17) 

f A - 

f o C A 

(15) 

r DA 

* r Do C A 

(18) 

REDUNDANTS 






6 a 

XX 

T 

f x M x 

(19) 

& xT 

a f T e_ 
x T 

(21) 

6 a 

XO 

o 

+ 

o 

3 

(20) 

6 xR 

" \ * r Dx ^D 

(22) 


1 ■ < B *o * * 8 *r - V (®3) 

ELEMENT FORCES AND STATICALLY DETERMINATE FORCES 


’x 1 + f o * 

(24) 


(25) 


DEFLECTIONS 

* ■ < f I 0 * + <4 D Fo * D to> f* < f I »T * V ' -DA (26 > 
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or 

A - » + V p + (F I D Po ♦ ^ * < F I e T + V - "L - X l \ < 2 ?> 


hereafter as statically determinate forces. The redundants are also unknown, of 
course, but the statically determinate forces resulting from unit values of the 
redundants are calculated, and these results are used in the continuity analysis. 
Calculating the statically determinate forces resulting from unit values of the 
external loads is also expedient. 

Each of the three sets of forces - external, redundant, and statically de- 
terminate - is further divided into two subsets called principal and subordin- 
ate forces. The subordinates are forces which can be expressed in terms of the 
principals by a preliminary calculation performed on the machine, after which the 
subordinate forces are eliminated from the problem. The principal forces are the 
forces that remain. The purpose of this elimination iB to conserve machine ca- 
pacity. 

The choice of subordinates should be such that they can be expressed easily 
in terms of their principals. For example, consider the shear panel of Figure 4. 
The forces on this panel form a self-contained system, and any three can be writ- 

«* a/b. The 

force is the principal, and Q gl , Q g2 , and Q g ^ are subordinates. By this device 
often half of the forces can be eliminated from the problem. 

The next step in the analysis, then, is to designate and number consecutively, 
beginning with one, each of the following six sets of forces: principal and sub- 

ordinate statically determinate forces, redundants, and external loads. Matrices 

of these forces are denoted respectively by the symbols Q, Q , X, X , P, and 0 . 

s s s 

Figure 3 illustrates a typical free body diagram with the forces numbered. On 
this diagram, only statically determinate forces are shown. Redundants and ex- 
ternal loads are shown on separate sheets to avoid confusion. Principal force 
numbers are enclosed in parentheses; subordinate force numbers are not. 

After the principal and subordinate forces are chosen, so-called subordinate 
degree of freedom vectors corresponding to equations of equilibrium are assigned, 
so that the subordinate forces can be calculated in terms of their principals. 


ten in terms of the fourth. Thus Q gl = a/b, Q g2 * and Q 
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These vectors are shown in the figure by dotted arrows with index numbers not en- 
closed in parentheses. Finally, principal degree of freedom vectors are assigned 
to permit the calculation of the principal statically determinate forces. The 
principal degree of freedom vectors are indicated by dotted arrows with index 
numbers enclosed in parentheses. In general, four sets of degree of freedom vec- 
tors are assigned as follows: principal degrees of freedom, and subordinate 

statically determinate, redundant, and external load degrees of freedom. 

The equilibrium equations can now be written, in matrix notation, in terms 
of the six sets of forces acting on the free bodies, as follows: 



The forces acting on the free bodies are contained in the post multiplier; 
the generalized components are contained in the preraultiplier . The significance 
of the partitions m^, etc., is given in detail in the table of notation. 

All of the generalized components are computed by equations 1, 2, 3# and 4. The 
null partitions in the generalized component matrix result from choosing subordin- 
ate forces in such a way that they always form small self-contained systems with 
their principals. 

Equation 28 is expanded as follows; 


m Q + m Q + p 
PP ps *x 


14 


0 


(29) 



in Q ■Hu Q n 
sp SB s 


p x x + p . x . 

sp ss 


0 

0 


P 0 + P 0 "0 

*0 v 8 

sp ss 


Q„ 


-mi m Q 
ss sp 


K 


w 

-P p JC 

_ *x 


ss sp 

-1 

ss sp 


“ p o V ^ 


The matrices IC^, K x , and K 0 are now defined according to equations 5, 6 , and 
7 of Table 1. 





Substituting these expressions into equation (29) gives 


MQ + P X + P j 6 - 0 , (30) 

x o r 

where the matrices M , P x , and P q are defined by equations 8 , 9> nnd 10 of Table 1. 
Equation ( 30 ) is the principal equilibrium equation. 

Notice that the matrices m , p , and p , appearing in equations 5, 6, 

88 x ss °ss 

and 7 t must be nonsingular. This nonsingularity is obtained by proper choice of 

subordinate degree of freedom vectors. As a matter of computing convenience, the 

choice of these vectors should be such that the matrices m , p , and p are 

88 x ss ss 

lower triangular, because in this event a very rapid computing program can be used 
to solve the equations. Such a choice is always easy to make, and it has the ad- 
ditional advantage that a lower triangular matrix with nonzero elements everywhere 
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on the diagonal is nonsingular, and veil -conditioned. 


Taking X ■ 1 (the unit matrix) and 0 «= 0 (the null matrix) in equation ( 30 ) 
leads to equation (ll) of Table 1, where is a matrix of statically determinate 
forces resulting from unit values of the redundants. 

Taking X = 0 and 0=1 leads to equation (12), where is a matrix of stati- 
cally determinate forces resulting from unit values of the external loads. 

Check degree of freedom vectors are assigned to various free bodies of the 
structure so that additional check equations are generated. Such equations pro- 
vide reliable verification of the calculations up to this stage. 

After the equilibrium problem is solved and checked, two additional opera- 
tions are performed, before the continuity of the structure is restored. First, 
all of the statically determinate forces, the redundants, and perhaps some of the 
external loads, are grouped into a single set of forces, called element forces, to 
facilitate calculating deflections. Secondly the statically determinate reactions 
are grouped into a separate matrix, to permit calculating the effect of support 
displacements . 

Element forces are defined in the following way; Consider any element of the 
structure which is capable of undergoing deformation, and therefore of contribut- 
ing to the deflection of the structure as a whole. Both internal forces and ex- 
ternal loads may act upon such an element, since the possibility of external loads 
acting between Joints is not excluded. Certain forces acting on the element are 
designated as element reactions . These element reactions may be internal forces 
or fictitious forces, but they must be chosen in such a way that they are capable 
of balancing the other forces applied to the element. The remaining internal 
forces are designated as element forces . After element forces for the entire 
structure are selected, they are numbered consecutively beginning with one. 

For each element force there is a corresponding element deformation . An ele- 
ment deformation is defined as the component of the displacement of an element 
force, in the direction of the element force, when the element reactions are un- 
displaced parallel to themselves. 
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Figure 5 shows a bending element, with element reactions (indicated thus? 
—»—*►), element forces (F^, F g , F^), and element deformations (e^ e 2 , e^). 

Other choices of element reactions, forces, and deformations are possible for 
such an element. 

The element deformations are given the same index numbers as the correspond- 
ing element forces; and a deformation is positive when it has the same direction 
as a positive value of the corresponding force. The sign convention for element 
forces is arbitrary, except that the choice of a sign convention which results in 
negative off-diagonal flexibility factors (defined later) is not advisable. 

Some of the element forces correspond to statically determinate forces; 
others correspond to redundante and a few may correspond to external loads. 
Therefore, the element forces can be written in terms of the statically determin- 
ate forces, the redundants, and the external loads, as follows? 

F-BQ + 1^X + H O 0, ( 31 ) 

where F is a matrix of element forces. 

If the element forces have been chosen in such a way that each one corres- 
ponds exactly to a statically determinate force, a redundant, or an external load, 

and such a choice should be made, then the matrices N. H , and H contain l's and 

’ x' o 

0's, and there will be no more than one 1 in any row or column. Such matrices are 
called extractors, because their only function is to extract information from 
other matrices. 

Setting X » I and fi = 0 in equation (31 ) yields equation (13) of Table 1, 
where f is a matrix of element forces resulting from unit values of the redun- 
dants. Setting X «* 0 and => I yields equation (l4), where f is a matrix of 
element forces resulting from unit values of the external loads. 

In the Maxwell -Mohr method, deflections are calculated by applying unit dum- 
my loads coinciding in position and direction with the desired deflections. In 
the present formulation the assumption is made that a unit external load is ap- 
plied to coincide with every such deflection. Therefore, a matrix f can be ex- 
tracted from f Q , as in equation (15), where f^ is a matrix of element forces 
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resulting from unit values of the dummy deflection loads, and is a suitable 
extractor matrix. 

Number the statically determinate reactions consecutively beginning with 1. 
Then the statically determinate reaction matrix can be extracted from the 
statically determinate force matrix as follows: 

«D " N r«' 

where N r is a suitable extractor. Setting X and 0 equal to I and 0 in turn leads 
to equations (l 6) and (17), where and r Do are matrices of the statically de- 
terminate reactions resulting from unit values of the redundants and external loads 
respectively. A matrix r of statically determinate reactions resulting from unit 
values of the dummy deflection loads is extracted from r_ as in equation (l8). 

The essentials of the derivation of equations (19) to (26), inclusive, have 
been given in reference 2. A feature of this derivation is that although it is 
based on the conservation of energy, it does not involve elastic strain energy, 
so that the deflection equations are immediately valid for arbitrary element de- 
formations, including deformations resulting from thermal gradients, plasticity, 
creep, etc. The derivation is also facilitated by the use of the notions of ele- 
ment reactions, forces, and deformations, as defined above. However, the equations 
have been generalized to include the effects of support displacements, the applica- 
tion of external loads between Joints, and the calculation of deflections at points 
between Joints. 

The symbol D appearing in these equations denotes the flexibility matrix. 

The elements of this matrix represent element deformations resulting from unit 
values of element forces. For example, the flexibility coefficients for the beam 
element of Figure 5 are as follows, if shear deformations are not considered: 

D L1 - L/AE, D 22 - L 3 /3EI, D 23 « D 32 - L 2 /2EI, - L/EI, 

where L, A, I, and E are the length, area, moment of Inertia and modulus of elas- 
ticity of the member. 

The matrix Dp Q contains element defo.. nations resulting from external loads 
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applied directly to the elements. If loads are applied only at Joints, then 
la null. Figure 6 shows the element of Figure with an intermediate load. 


The following elements of the D„ matrix can he derived hy elementary methods i 

J? o 

D„ - a cos a/AE, D„ » a 2 (L-a/3) sin a/2EI, D - a 2 sin a/2EI . 

Fo lJ Fo 2j i0 3J 


The matrix D contains displacements of dummy deflection loads acting direct- 
AF 

ly upon the element, resulting from unit values of the element forces, when the 

element reactions are not displaced parallel to themselves . Figure 7 shows the 

element of Figure 5 with an intermediate dumay deflection load. The elements of 

D a _ are bb follows! 

AF 


AF 


b cos 


J1 


l /AE, 


AF 


J2 


(L-b/3) Bin^/2EI, 


AF 


- b" 


J3 


sinjS/2EI 


The matrix contains displacements of dummy deflection loads acting direct- 
ly upon the element, resulting from unit external loads acting directly upon the 
element, when the element reactions are not displaced. Figure 8 shows a bending 
element subjected to intermediate external and deflection loads. The correspond- 
ing element of is as follows: 

2 

D Ao “ |a cos a cos p > + Bin a 6ln 

i J 

or 

2 

D Ao * Ia COS ° cos f + sin a sin ^ if b)>a . 

i J 

The matrix e^, contains element deformations resulting from heating, plastic- 
ity, creep, etc. For example, suppose that the tensile element of Figure 9 (») 4 
has been assigned the ith element force, as shown. In (b) the temperature of the 
element is increased an amount AT in the Jth analysis condition. The thermal de- 
formation is then e., «* a LAT, where a Is the coefficient of expansion. The 

X iJ 

matrix elements e can also represent bending thermal deformations of bars heat- 

U 

ed unequally on the two sides, or any other kind of a thermal deformation. When 

the e_ represent plastic or creep deformations, they either must be known, as 
X i 
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they could be In a statically determinate structure, or they must have been com- 
puted in a previous cycle of some kind of iterative process. 

The matrix e^ contains displacements of the dummy deflection loads acting 

directly upon the element, resulting from heating, etc., when the element reactions 

are not displaced parallel to themselves. Figure 10 shows the element of Figure 

9 with an intermediate dummy deflection load. The intermediate thermal deformation 

is e „ = a a AT. 

m iJ 

The matrices and contain displacements of the statically determinate 
and redundant reactions, respectively. The elements of these matrices are posi- 
tive when the corresponding support displacements have the same sense as positive 
values of the reactions acting upon the structure. 

Equation (27) provides an alternate, more accurate, but somewhat more cumber- 
some means of calculating deflections. In this equation, F^, and R DA are ma- 
trices containing element forces, redundants, and statically determinate reactions, 
respectively, in the uncut structure resulting from unit values of the dummy de- 
flection loads. The equation can be shown to be mathematically identical to equa- 
tion (2 6). 


COMPUTER PROGRAMS 

The calculations are performed on an IBM 709 computer. The only " 709 " pro- 
gram written specifically for the Maxwell -Mohr method is called "Matrix Generation" 
This program accepts, as input, coordinates and directions numbers which define the 
degree of freedom and force vectors appearing on the free body diagrams. The di- 
rection numbers have previously been computed from the coordinates by an auxiliary 
program. Thus, the only numerical input prepared by the analyst for this phase is 
a table of coordinates. The program then generates the elements of the matrices 


“pp* m ps* p x 1 p x * p o * p o * “sp* m ss' P x 

pp ps pp ps * sp 

of equations (l) to (4) of Table 1. 


> P x > P Q , and p Q by means 
SS sp SB 


All the rest of the calculations, as required by equations ( 5 ) to ( 26 ), are 
performed with the aid of a general purpose interpretive routine called the "Tape 
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Matrix Compiler". This routine essentially permits the analyst to write his own 
programs for matrix operations. Matrices of member flexibilities, loads, thermal 
deformations, and support displacements, and certain extractor matrices, are in- 
put. The machine outputs the unknown forces and deflections of the structure. 

The compiler is also used to perform additional operations not covered by 
equations (5) to (26). These auxiliary operations can include transforming the 
redundants to improve conditioning, and the modification of member flexibilities, 
including the complete removal of members. 

The Joining of structures to form larger structures iB accomplished by the 
basic program, comprising equations (l) to (26). 

A program under development, called the "Structure Cutter", permits the 
machine to select its own redundants optimized to yield well -conditioned equa- 
tions. The capabilities of the Structure Cutter are briefly discussed in a later 
paragraph. 
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ANALYSIS PROCEDURES 


IDEALIZINO THE STRUCTURE 

The actual structure is replaced by an idealized discrete structure 
consisting usually of bars and panels. In general the bars can carry tension, 
torque, two components of bending moment, and two components of shear. The 
panels can carry shear and biaxial tension. In the most generally useful 
idealization, bars are considered straight between joints, and panels carry 
only shear. However panels are permitted to be warped » This allowance for 
panel warping Improves the accuracy of the analysis, because joints of the 
idealization can lie on the true contour of the actual structure. Furthermore, 
warping simplifies the input, because there are few if any derived coordinates. 

The meaning of the term "derived coordinates" is explained later. 

Panels should be rectangular if possible, trapezoidal if not rectangular, 
or at least nearly trapezoidal. Panels that almost come to a point should be 
avoided. Triangular panels should probably be removed, leaving a triangular 
framework of bars. 

A problem of structural idealization concerns the question of the attachment 
of shear panels to bars. Two methods of attachment are considered. In the first 
method, panels are attached to bars at the midpoints of panel edges, as shown at 
"A" of Figure 12. In the second method, the attachment is continuous, as shown 
at "B", and the assumption is made that load in the adjacent bars varies linearly 
between Joints. 

Figure 11 shows a set of skin-stringer panels, rigidly supported at infinity. 
The panels have symmetry about the X-axis, the stringers are equally spaced and 
have constant area, all the stringers are equally stiff, and the sheet thickness 
is constant. Transversly the panels are assumed to be stiffened by a continuum 
of infinitely rigid bars. Axial loads are applied to the #3 stringers at X ** 0. 

The exact solution of the stringer loads and panel shear flows in the struc- 
ture was obtained. The structure yati also analyzed by the Maxwell-Mohr method, 
for the idealization shown in Figure 12. At X «= 80, conditions are essentially 
the same as they are at infinity. 
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Two digital solutions were obtained. In the first solution, panels were 
assumed to be attached to bars at panel mid-points only. Under this assumption, 
the load in a bar is constant, but can jump abruptly at Joints and panel mid- 
points. The flexibility matrix corresponding to this assumption is diagonal. 

In the second solution the load in the bar is assumed to vary linearly 
between Joints. The flexibility matrix in this case is not diagonal. 

The comparison of the three solutions for stringer loads is given in Table 
2. The results for methods 1 and 2 are followed by the percent errors in pa- 
rentheses. The comparison for shear flows Is given in Figure 12. 


Table 2 Stringer Loads 


Strlnger\ X 
Number \ 

0 


20 

40 

60 80 


Exact 

0 


.1470 

.1889 

.1987 .1996 

1 

Method 1 

0 

(0) 

.1446 (-2 j>) 

.1869 (-1$) 

.1970 (-1 $>) .1988 (-0#) 


Method 2 

0 

(o) 

.1344 (-9$) 

.1885 (-0$) 

.1986 (-056) .1998 ( 0 i) 


Exact 

0 


.2015 

.2030 

.2007 .2002 

2 

Method 1 

0 

(0) 

.1925 (-45Q 

.2020 (-056) 

.2008 ( 0j>) .2004 ( 056) 


Method 2 

0 

(0) 

.1864 (-8$) 

.2060 ( 1$) 

.2005 (-056) .2001 (-0?6) 


Exact 

1.0000 

.3031 

.2152 

.2029 .2005 

3 

Method 1 

1.0000 (o) 

.3258 (7#) 

.2221 ( 3 1=) 

.2044 (l/o) .2017 (l$) 


Method 2 

1.0000 (o) 

.3584 (18 io) 

.2110 (-2$) 

.2018 (-1$) .2004 (-0$) 


The comparisons show that the "panel mid-point method" gives greatest accuracy. 
The fact should be noted however that this method gives somewhat less accuracy than 
the second method for the deflection of a cantilever thin web beam, idealized as 
shown in Figure l4. Here the accuracy of the deflection computed by the first 
method depends on the number of bays and is satisfactory for four bays. Both 
methods give correct cap loads and shear flows for any number of bays. 

Since the "panel mid-point method" is the simplest, and seems to be the 
most accurate, at least for stresses. It appears to be preferable to the second 
method. 
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A distinction is made between "defining" and "derived" coordinates. This 
distinction is demonstrated in Figure 15, which shows a pin-jointed truss lying 
in the X-Y plane. Member AC is assumed straight. The geometry of the truss 
therefore may be considered to be defined by the X and Y coordinates of points 
A, B, and C, and the X coordinate of D. The coordinate Y^ can be derived from 
Xp on the assumption that AC is straight. The coordinate Y D is therefore a 
derived coordinate, and the others are defining coordinates . Defining coordi- 
nates should be input with an accuracy of about six decimal places to avoid 
contradictions between them and the assumptions upon which they are derived, 
within the machine. Because of this accuracy requirement, derived coordinates 
should be avoided. 

A warped shear panel cannot be in equilibrium under the action of shear 
forces alone, as Figure 16 demonstrates. The shear forces shown in the plan 
view all have downward components in the edge view. The panel can be put into 
equilibrium with the addition of two forces at opposite corners, as shown in 
the perspective view of Figure 17 . This figure also shows principal and sub- 
ordinate force numbers, and subordinate degree of freedom vectors, which can 
be assigned to permit the machine to calculate the subordinate forces in terms 
of their principals. The warping forces are approximately normal to the panel. 

The reactions to the warping forces are assumed to act on joints. 

Many structures contain warped panels which cannot be flattened in the 
idealization without seriously compromising the accuracy of the solution. 
Furthermore, the flattening process is usually more trouble than accounting for 
the warping. 

CUTTING THE STRUCTURE 

Box structures, like wings, composed of bars in tension and panels in shear, 
tend to be better conditioned, because they are stiffer, than fuselage-type 
structures which contain flexible rings. For structures which are inherently well 
conditioned, and yet which may offer cutting difficulties because of unusual 
features, the "building method" is a useful procedure. 

In the building method, a unit of the structure known to be statically 
determinate is selected, and the structure is built from this unit by adding 
other statically determinate units. The members which are omitted in the process 
are the redundantc. 
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Figure 18 (a) shows the uncut structure, (b) and (c) si stages in the 
building method and (d) shows the final cut structure. Two panels and a re- 
action are redundant. In the process, the use of "temporary reactions" may 
be expedient. These reactions can be replaced by the actual reactions at the 
completion of the process. 

The following expression is convenient for checking the degree of redundance 
of a structure composed of shear panels and axially loaded bars? 

n«*b+p + r- 2Jg - where b » the number of uncut bars, 

p m the number of uncut panels, 
r ■ the number of reactions, 

Jg * the number of two constraint Joints, 

J^ •» the number of three constraint joints, 

For a statically determinate structure, n « 0. The expression, with n = 0, is 
a necessary but not a sufficient condition for static determinacy. For the 
structure of Figure 18, n=28 + lk + 6- 2x0-3xl6»0. 

DIAGRAMS 

The following diagrams are utilized: (l) a general view of the Idealized 

structure with the joints numbered consecutively beginning with one, (2) a set 
of free body diagrams, and ( 3 ) diagrams showing the element forces. 

The free body diagrams have been described in the section on matrix formu- 
lation, and Figure 3 shows a typical diagram for statically determinate forces. 

The only feature of these diagrams not already mentioned are the free body numbers, 
shown enclosed in squares in Figure 3» The machine uses these numbers to associate 
forces with their corresponding degrees of freedom. 

The element force diagrams show element reactions and element forces, the 
latter being numbered consecutively beginning with one. The statically determi- 
nate forces and redundants should he chosen so that, each element force is identi- 
cal with either a statically determinate force or a redundant, so that the elements 
of the If and matrices consist only of l's and 0‘s. 

LOAD SHEETS 

Data is input on three different formats, as follows: the coordinate table, 
the vector description tables, and the matrix load sheet. 
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The coordinate table is a list of Joint numbers with their associated X, 

Y, and Z coordinates. With the aid of an auxiliary program, the machine com- 
putes a table of direction cosines of vectors defined by point pairs of the 
coordinate tables. The point pairs are specified by the analyst on a separate 
load sheet. The auxiliary program can also compute the direction cosines of a 
vector defined as the cross product of two other vectors each in turn defined 
by point pairs designated by the analyst. The vectors for which direction 
consines are calculated include most, or all, of the vectors which appear in 
the analysis. Direction numbers of additional vectors can be hand input if 
necessary. The machine sorts the computed direction cosines according to the 
defining points, and assigns each set of X, Y, and Z direction cosines a serial 
number. 

The vector description tables are of two types. On the type 1 table the 
following information is input for each vector: the vector serial number; the 
type, whether angular or linear; the sign; the number of the free body upon 
which the vector acts; the number of a point on the line of action of the vector; 
and the serial number of the direction of the vector. Each vector is listed 
only once in the type 1 load sheets. However most of the force vectors appear 
more than once on the free bodies, and an entry must be made each time a vector 
appears. These additional entries are made on the type 2 tables which have 
provision only for vector serial numbers, signs, and free body numbers. The 
type 1 and type 2 tables are filled out for the four kinds of degree of freedom 
vectors, and the six kinds of force vectors mentioned previously. 

The matrix load sheets contain spaces for the matrix elements, and for 
the row and column numbers corresponding to each element. The matrices If, H^, 
H o , C a , D, D f0 , D ao , e T , e^, Ap, A x and are input on these sheets. Occasion- 
ally some elements of the K q matrix also are hand input. 

Ordinarily only the matrices H, H^, H q , C a , D and <t> are required, and of 
these matrices H, H v , H and C should contain only l's and 0's. Thus the only 
formats which contain numerical input are the coordinate table, the flexibility 
matrix D, and the load matrix Therefore a problem which. jhas been set up for 
a given sot of coordinates, flexibilities, and external load can be solved for 
new coordinates, flexibilities, and loads by inputting only three tables. These 


2 6 



tables represent the minimum possible input for the problem. Therefore a given 
set-up, say for a fuselage section, can be used many times for a variety of 
fuselage analyses, and the set-up essentially becomes, in itself, a general 
program for fuselage problems. 

CHECKS ON THE OUTPUT 

The equilibrium checks, made by writing extra equations of equilibrium, 
have been mentioned. Two other important types of checks are the simultaneous 
equation checks and the symmetry checks. Simultaneous equation checks are made 
on the solutions of both the equilibrium and the continuity equations by sub- 
stituting the results into the original equations. A symmetry check is made on 
5 , which must be symmetric by Maxwell's law. A similar check is made on the 

deflection matrix A, for rows and columns which correspond to identical unit 
deflection loads and external loads. 

IMPROVING THE CONDITIONING 

Naturally every effort should be made at the beginning to secure well- 
conditioned equations. The familiar rule la that redundants should be chosen 
which are small compared to other forces in the structure. The rule can alBC 
be stated as follows: in the cutting process the structure should lose as 

little stiffness as possible. For example, a good choice of redundants for a 
fuselage frame is the insertion of three hinges . A complete cut at one point 
leaves the frame very flexible. 

A second device Is to break the structure into statically indeterminate 
substructures.. The substructures are then cut and analyzed, after which they 
are joined to form the original structure, as discussed in a later paragraph. 

At each stage of this process the redundants are relatively few in number, and 
generally well conditioned. 

A third device is the utilization of orthogonal solutions derived from the 
theory of elasticity for geometrically regular bodies which resemble the struc- 
ture at hand. This process has been thoroughly discussed by Argyris. 

JOINING SUBSTRUCTURES 

In this process the structure is broken, by cutting redundants, into sub- 
structures, which remain joined together by other forces which can be computed 
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from statics. Thus the cut structure can he regarded as n statically determinate 
structure consisting of statically indeterminate elements. Figure 19 shows a 
DC-8 wing-pylon intersection which has teen broken into two substructures by this 
method. The figure shows statically determinate forces only. The other joining 
forces, which are redundants, are shown on a separate sheet. Figures 2a and 2b 
show details of the idealized substructures. 

After the structure has been cut into substructures, each of the substructures 
is also cut and analyzed in detail, for unit values of the external loads, which 
include the Joining redundants. In particular the deflections of the substruc- 
tures, at points where they have been cut apart, are calculated. The analysis of 
each substructure utilizes the basic program and the equations of Table 1. 

After the substructures are analyzed, they are joined to form the original 
structure by another application of the basic program and equations. In this 
process free body diagrams, like figure 19, are drawn. Element force diagrams 
are also prepared. Element reactions for the substructures, considered as 
elements of the original structure, must be identical with the statically de- 
terminate reactions that were utilized in the detail analysis of the substructures. 
This requii'ement is necessary because the elements of the flexibility matrices 
D, D._, and D A ^ are extracted from the deflection matrices A, calculated for 

each of the substructures. The extraction is accomplished with the aid of ex- 
tractor matrices consisting of l's and 0's and the tape matrix compiler. 

DISC 0NN3CTIHG AIID FLEXIBILITY MODIFICATION 

The technique discussed by Argyris , Michielsen and Dijk [l3] > and 
Best Ql 4J , for modifying flexibilities with the aid of arbitrary element defor- 
mations after the redundants have been computed , has two important applications . 
First the effect of changing the sizes of a few members upon the stress distri- 
bution can be determined with a minimum amount of calculation. However the method 
becomes Inefficient when the number of elements to be modified becomes equal to 
or greater than the number of redundants. In this case a new flexibility matrix 
should be input. Second, the notion of filling in cut-outs, like fuselage doors, 
and later removing them, is important, because the process of cutting the struc- 
ture is greatly simplified when cut-outs are not present, and the equations are 
likely to he better conditioned. However, more machine capacity is required. 
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Members can also be removed by making them more flexible, say on the order 
of 1,000,000 times, than other members of the structure. This approach only 
works when the forces being reduced to zero are redundants. Otherwise the 
continuity equations tend to be linearly dependent. 

THE STRUCTURE CUTTER 

A method has been devised for having the machine cut the structure. In 
this approach no distinction is made between statically determinate and re- 
dundant forces when the problem is set up. The number of unknowns in the equi- 
librium equations generated by the machine then exceeds the number of equations. 

By a process of selecting columns of the rectangular matrix of coefficients of 
unknowns in these equations, the machine chooses a well-conditioned square matrix. 
The unknowns which correspond to the columns of this matrix are the statically 
determinate forces, and the remaining unknowns are the redundants. The choice 
of columns is influenced by weighting factors which reflect the stiffness of 
the members of the structure. 

Figure 20a shows a statically indeterminate structure. Figure 20b shows 
the same structure as it was cut 'by the machine. 

SIMPLIFIED IHPUI 

A new program. called the "Redundant Force Method" is being developed. This 
program is basically the same as the method described previously, but the new 
method incorporates a number of improvements which eliminate the need for pre- 
paring free body diagrams, and reduce the input to a minimum. In effect the 
machine automatically cuts the structure (utilizing the "Structure Cutter"), 
breaks the statically determinate structure into free bodies, writes and solves 
the equations of equilibrium, and writes and solves the equations of continuity, 

A certain penalty in additional machine time is involved, however the new program 
is expected to be especially useful in the rapid solution of preliminary design 
problems for which a rough idealization is satisfactory, and which cannot be 
solved without a large error by elementary methods. 

NOHLIIIEAR PROBLEMS 

Although tnls subject is beyond the scope of the present paper, some mention 
should be mode of the applications to the nonlinear problems involved in calculating 
the effects of plasticity and creep upon the behavior of the structure. The approach 
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to these problems has been through the use of various step-by-step, or iterative, 
procedures. In all such procedures the question of convergence is of primary 
importance, because the rate of convergence can be fast or slow, or the process 
can be divergent. Rapid convergence iB necessary, because a large amount of 
calculation per cycle is required even for a structure of moderate size. 

A method of calculating stresses and deflections in the presence of plasticity 
is given in reference 6. The method utilizes the rapidly convergent Newton-Raphson 
procedure for solving nonlinear simultaneous equations. Agreement with test results 
is demonstrated. Reference 7 presents an approach based on the use of fictitious 
loads which appears to require a minimum amount of computation per cycle. 

A step-by-step application of the Maxwell-Mohr analysis to the creep problem 
is under development. This work is expected to provide a means of computing the 
history of stress and deflection of a statically indeterminate structure subjected 
to time dependent load and thermal inputs. 
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COMPARISON WITH TEST RESULTS 


Comparisons between analysis and test results obtained at the NASA and 
during the DC -8 static test have been made. The NASA comparisons were ac- 
complished in the period from June 1956 to September 1957* In a ll the numeri- 
cal analysis, the midpoint idealization for shear panels was used. 

The comparison for axial stresses measured in the cantilever circular 
cylinder of Figure 1 has been mentioned. Figure 21 shows the analytical and 
test results for frame bending moments and skin shear flows in the same cylinder. 
The results of the Maxvell-Mohr analysis are in very close agreement also with 
results obtained by the method of Hoff , as reported in reference 9 . 

Figures 22, 23, 24, and 25 show comparisons for a swept box tested at the 
NASA, and reported in reference 10 . The box was of rectangular section and 

had a total of 32 stringers. In the figures the heavy solid lines indicate 
idealized stringers and bulkheads, while the dotted lines indicate bars obtained 
by lumping skin in the chordwlse direction. The analysis would not yield sat- 
isfactory approximations for shear flows in the covers until these bars were 
inserted. Poisson's ratio was accounted for in the triangular area at the root. 
In the bending test, the characteristic peaking of axial stress at the rear spar 
is correctly predicted, as is the reversal of shear flow in the front spar web. 

Figures 26 and 27 show comparisons for cylinders with cutouts subjected to 
bending and torque respectively. The tests are described in references 11 and 
12. As the figures show, more idealized stringers were inserted in the upper 
aide than in the lower, because the cutout at the top perturbs the stress field, 
and requires finer lumping. Frame flexibility was taken into account. The 
resulting agreement is excellent. However there is one shear panel at the 
corner of the cutout which, in the bending case, does not have approximately a 
uniform, shear flow, as assumed. At one edge of this panel the shear flow, not 
shown in the figure, is considerably higher than the value at the panel center. 
The only way to cover this concentration without going to a finer lumping is 
with an empirical factor. 

Figure 20 shoves a comparison of measured and calculated stresses for a 
station in the root region of the DC-0 wing. The analysis which yielded the 
calculated results is discussed in a later section. 



APPLICATIONS 


The method has been extensively applied in the analysis of Jet transport 
components; missile parts, including fins and body components; and a supersonic 
low aspect ratio wing. Many of these analyses included calculations of thermal 
stress and deflection. 

The wing-fuselage intersection was one of the primary problems in the stress 
analysis of the DC-8. The stress distribution was complicated by the existence 
of wing sweep, an auxiliary spar, landing gear cutouts in the lower part of the 
fuselage behind the wing, a keel beam running along the fuselage centerline 
below the floor, and other details. . The problem was approached by first making 
an analysis of the entire region, including a fairly detailed representation of 
the fuselage, and a simplified idealization of the wing. From the results of 
this analysis, reaction forces between wing and fuselage were determined. A 
detailed wing root analysis was then made, In which these reaction forces were 
applied . 


Figure 29 is a diagram of the idealised structure used in the detailed wing 
analysis, showing the three spar conatructissss# with the auxiliary spar which 
supports the main landing gear. The idealization had the correct sweep; dihedr al ; 
incidence and taper, both in plan-form and in thickness; and the airfoil sections 
were accurate. However, twist was removed to flatten skin panels. There were 
113 redundants and 300 element forces. The first complete calculation based on 
this idealization was finished in March 1956. Had the Job been done a little 
later, panel warping and twist would have been considered. 

The idealized structure for the tail- fuselage intersection is shown in 
figure 30. The idealization Included a portion of the vertical tail, and a 
stub of the all-movable horizontal surface. Some of the sections were stiffened 
by frames like the one shown in section A-A; others had partial bulkheads. The 
Joints of the idealized structure lay on the true contour, and panel warping was 
accounted for. The foreward and aft parts of the structure were analyzed sepa- 
rately and then Joined at section A-A. The first complete calculation was 
finished in September 1957. 

Deflection Influence coefficients calculated for both the wing and the 
fuselage tail section wore used in flutter analysis. 
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An analysis of the fuselage nose section, including the cockpit enclosure, 
was performed. The problem was complicated by the presence of cabin pressure, 
and the fact that the pressure envelope was irregular because of the existence 
of the unpressurized nose-wheel well below the floor. The members of the cock- 
pit canopy also caused added difficulties, because some of them were designed 
to carry tension, bending moments about two axes, and torque, and they were so 
analyzed. The structure was analyzed in two separate sections, which were then 
Joined. After joining, the technique of virtual disconnecting loads was employed 
to calculate the effect of door cutouts. 

Figure 2a and 2b show the idealized structure for the Conway outboard pylon. 
The structure incorporates a bottoming strut, shown in figure 2a. The bottoming 
of this strut, after a certain amount of load has been applied, changes the stress 
distribution, and causes a nonlinearity, which wa s taken into account. 

Figure 31 shows the structure of the JI-4 ejector-reverser. The structure is 
irregular; has large cutouts for the reversing buckets; incorporates members sub- 
jected to tension, bending about two axes, and torque; and is subjected to large 
thermal gradients. The JT-3 and Conway ejectors are similar. Results from the 
JT-3 analysis became available within a period of two, months. The same set-up 
was then utilized in the analysis of the JT-4 and Conway ejectors, which have 
different sizes, shapes, and stiffnesses. The Conway ejector analysis was com- 
pleted in final form ready for submission to the FAA in one month's time. Spring 
constants for the ejectors were calculated and shown in proof test to be correct 
within the experimental error. 

Numerous applications to low aspect ratio wing and missile structures have 
been made, but these projects are classified and cannot be discussed. However 
the foregoing applications and experimental verifications have demonstrated that 
the matrix equations and the computer program are sufficiently general to deal 
with any linear discrete structure. Missile and supersonic airplane structures 
are no exceptions. Thus the low aspect ratio multi-spar wing- fuselage structure 
of figure 32 can be analyzed, with all the detail shown and more, with Joints on 
the true contour, for load and thermal stress. Deflections, and a deflection 
influence matrix useful in flutter analysis also can be output. 
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CONCLUSION 


A procedure for structural analysis, comprising a matrix formulation of 
the equilibrium and Maxwell-Mohr continuity equations, and an associated digital 
computer program, has been developed. This procedure is applicable, in its 
basic form, to any linear discrete structure. The method has been fully veri- 
fied by comparison with test results, both in the laboratory and in proof test, 
and it has been shown to be a practical analysis tool in numerous applications. 

Procedures of this kind, several of which have appeared in the last few 
years, represent a break-through in the art of stress analysis. These methods 
permit the practical calculation of stresses in complicated shell structures 
in rigorous accord with basic physical principals. This rigor is necessary, 
because approximate methods widely used in the past can be in error by large 
amounts. These errors are alleviated somewhat by stress redistribution above 
the yield, but below the yield they represent stress concentrations which cause 
premature fatigue failures. Above the yield premature static failures can occur 
in spite of the redistribution. 

In the past, serious consequences of these errors have been avoided by 
extensive testing. Some testing will always be necessary, but it is expensive, 
even compared to the cost of operating a large digital computer. In the future, 
testing expense will increase as airframes become larger, and the additional 
complication of thermal gradients is introduced. Therefore the need for rigorous 
methods is Increasing. 

Douglas Aircraft Company, Inc., 

Santa Monica, Calif., September 17, 1959. 
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Axial stress in a cylindrical Shell 
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Fig. 2a Idealized pylon structure 
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Fig. 2b Idealized wing structure 
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PLAN VIEW 
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Fig. 1 6 Warped shear panel 





Fig. 17 Equilibrium of a warped shear panel 
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Fig. 19 Free body diagram for joining structure 







Fig. 20a Statically Indeterminate structure 
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Fig, 20b Machine cut structure 
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NA.CA swept box 
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